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We give a pedestrian derivation of a formula of Biittiker et. al. (BPT) relating the adiabatically 
pumped current to the S matrix and its (time) derivatives. We relate the charge in BPT to Berry's 
phase and the corresponding Brouwer pumping formula to curvature. As applications we derive 
explicit formulas for the joint probability density of pumping and conductance when the S matrix 
is uniformly distributed; and derive a new formula that describes hard pumping when the S matrix 
is periodic in the driving parameters. 
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Brouwer ||, and Aleiner et. al. jij, building on re- 
sults of Biittiker, Pretre and Thomas (BPT) pointed 
out that adiabatic scattering theory leads to a geomet- 
ric description of charge transport in mesoscopic quan- 
tum pumps. Some of these works, and certainly our own 
work, was motivated by experimental results of Switkes 
et. al. j6| on such pumps. 

In this article we examine the formula of BPT ||, 
which relates adiabatic charge transport to the S matrix 
and its (time) derivatives, in the special case of single- 
channel scattering. We show that the formula admits a 
simple interpretation in terms of three basic processes at 
the Fermi energy. Two of these are dissipative and non- 
quantized. The third integrates to zero for any cyclic 
variation in the system. 

Next, we describe the geometric significance of BPT 
and relate it to Berry's phase JE|. It follows that the 
pumping formula of Brouwer P] can be interpreted as 
curvature and is formally identical to the adiabatic cur- 
vature ||] . In spite of the interesting geometry the topo- 
logical aspects of pumping are trivial. In particular, we 
prove that all Chern numbers associated to the Brouwer 
formula are identically zero. 

We proceed with two applications. First we give an 
elementary and explicit derivation of the joint probabil- 
ity density for pumping and conductance. This problem 
was studied in ||]. Brouwer 's results go beyond ours as 
he also calculates the tails of the distributions and we 
don't. On the other hand, parts of his results are nu- 
merical, and they are certainly not elementary. Finally, 
we calculate, for the first time, the asymptotics of hard 
pumping for S matrices that depend periodically on two 
parameters. If the system traverses a circle of radius R in 
parameter space, with R large, then the amount of charge 
transported is order VR, multiplied by a quasi-periodic 
(oscillatory) function of R leading to ergodic behavior. 

We shall use units where e = m = H = 1, so the elec- 
tron charge is —1 and the quantum of conductance is 
e 2 /h — The mutual Coulombic interaction of the 
electrons is disregarded. 



The BPT formula: Consider a scatterer connected 
to leads that terminate at electron reservoirs. All the 
reservoirs are initially at the same chemical potential and 
at zero temperature. The scatterer is described by its 
(on-shell) S matrix, which, in the case of n channels is 
an n x 7i matrix parameterized by the energy and other 
parameters associated with the adiabatic driving of the 
system (e.g. gate voltages and magnetic fields). 

The BPT formula § says that the charge dqi enter- 
ing the scatterer from the £-th lead due to an adiabatic 
variation of S is 

dq t = ^'Ti(Q i dSS^, (1) 

where Qi is a projection on the channels in the €-th lead, 
and the S matrix is evaluated at the Fermi energy. In 
the special case of two leads, each lead carrying a single 
channel, 

where r, (r') and t, {t') are the reflection and transmis- 
sion coefficients from the left (right) and Qi projects on 
the left lead. In this case Eq. (Q) , for the charge entering 
through the left lead, reduces to 

2irdqi = i{fdr + t'dt'). (3) 

We shall present an elementary derivation of (^) . 

Derivation: Every unitary 2x2 matrix can be ex- 
pressed in the form: 

_ i7 / oos(0)e fa isin(0)e-*\ 

\ism(6)e^ cos{8)e- ia J ' 1 ' 

where < a, <f> < 2n, < 7 < n and < 8 < n/2. In 
terms of these parameters, Eq. (||) reads 

2irdqi = - cos 2 (6) da + sin 2 (0)d<f> - dry. (5) 

The basic strategy of our derivation of Eq. (||) is to find 
processes that vary each of the parameters in turn, and 
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keep track of how much current is generated by each pro- 
cess. An underlying assumption is that current depends 
only on S(kp) and S(k F ), so that processes that give rise 
to the same change in the S matrix also give rise to the 
same current. Because we do not prove this assertion, 
our derivation cannot be considered a complete proof. 



2wdq = sin 2 '(6)d(j>. 



(8) 




FIG. 1. Moving the scatterer changes a — > a + 2k F dL 




FIG. 2. Applying a vector potential changes 
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We understand the four parameters as follows. (See 
figures.) The parameter a is associated with transla- 
tions: Translating the scatterer a distance dL — da/2k F 
to the right multiplies r, (r') by e lda , (e~ lda ), and leaves 
t and t' unchanged. The parameter <j> is associated with 
a vector potential A near the scatterer. This induces a 
phase shift d<f> = — f A across the scatterer, and mul- 
tiplies t, It') by e ld ' t ', (e~ 4<i ^), while leaving r and r' 
unchanged. The parameter 9 determines the conduc- 
tance of the system: g — \t\ 2 /2tt = sin 2 (#)/27r. Finally, 
7 = (— i/2) logdet S is related, by Krein's spectral shift 
0, to the number of electrons trapped in the scatterer. 
As a consequence, for any closed path in the space of 
Hamiltonians § d^ = 0. 

To determine the effect of changing a we imagine a pro- 
cess that changes a and leaves the other parameters fixed, 
namely translating the scatterer a distance dL = da/2kF 
to the right. The scatterer passes through a fraction \t\ 2 
of the k F dL/ir — da/2n electrons that occupy the region 



A current / then dissipates energy at the reservoirs at a 
rate 



P = 2wI 2 /\t(k F )\ 1 



(9) 



To understand the effect of changing 9 and 7, we first 
suppose our scatterer is right-left symmetric, so r — r' 
and t = t' . Then changes in 9 and 7 would draw equal 
amounts of charge to the scatterer from the left and right 
leads. The charge that accumulates on the scatterer is 
given by Krein's spectral shift |7l The charge coming 
from the left is half this, namely ||]: 
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2ndq — d logdet S = —dj. 

4tt 
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Since every S matrix can be obtained by translating and 
adding a vector potential to a right-left symmetric scat- 
terer, the formula ([l(]) applies to all possible iS-matrices. 
The formula (|Io| ) applies to all possible scatterers, sym- 
metric or not. 

Combining (§), (§) and ©, gives BPT, Eq. (§). 

The effect of changing 7 integrates to zero on a closed 
loop. Changing 9 does not give any transport at all. 
Thus, only changes in a and <fi contribute to the net trans- 
port of a quantum pump. These are dissipative processes, 
with the rate of energy dissipation P (in the reservoirs) 
that is bounded from below by 27r/ 2 . This is contrary to 
the assertions of who claimed that the charge trans- 
port is the sum of a quantized non-dissipative term and 
a dissipative term that is not quantized. 

Geometrical interpretation: A = —2irdq is the 1- 
form (vector potential) associated with Berry's phase [|| . 

If we define the unit spinor \ip) = I ^ ) then 
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i(ip\dip}- 



(11) 



of size dL, and pushes the remaining |r| 2 da/27r electrons T h e set of all spinors \ip) is a 3-sphere (since |r| 2 + \t'\ 2 = 



forward. Thus 



2-ndq = - cos 2 (9) da. 



(6) 



This result can be obtained less heuristically, || , by work- 
ing in the reference frame of the moving scatterer and 
integrating the contribution of each wave number from 
to k F - From this one also sees || that the rate of dissi- 
pation at the reservoirs, P, is quadratic in the current I, 
with a coefficient that depends on the dispersion relation. 
If the dispersion relation is quadratic, then 

P = 27rI 2 /\r(k F )\ 2 . (7) 

To change (j), we vary the vector potential. This in- 
duces an EMF of strength — J A — <j>. The current is sim- 
ply the voltage times the Landauer conductance \t\ 2 /2tt 
H . Integrating over time gives 



1), while the set of all ratios r/t' is the projective space 
CP 1 = C + {00} ~ S* 2 . The natural map between them, 
namely \tp) — > r/t', is called the Hopf fibration, and A is 
called the "global angular form" of this fibration. 

To compute the charge transported by a closed cycle C 
in parameter space, we can either integrate the 1-form A 
around C, or (by Stokes' theorem) integrate the exterior 
derivative (curl) f2 = dA over a disk D whose boundary 
is C. Q is the curvature 2-form of Brouwer [Hi: 



-i(dip\dip) = —%■ 



dz A dz 

1 + M 2 ) 1 



(12) 



where z — r/t' . The expression —i(dip\dip) is formally 
identical to the adiabatic (Berry's) curvature that ap- 
pears also in the quantum Hall effect [Q. 
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In the last expression one sees that the curvature sees 
only the ratio z = r/t' , and not r and t' separately. The 
curvature fl is the [7(2)-invariant area form on CP 1 , and 
its integral over all of CP 1 is 2n. Q is also the curvature 
of the Hopf fibration. 

In the study of non-dissipative quantum transport, 
Chern numbers play a role. These are topological invari- 
ants that equal the integral of the curvature over closed 
surfaces in parameter space. In the context of adiabatic 
scattering, however, all Chern numbers are zero. The 
vector bundle over parameter space is topologically triv- 
ial, and the vector (r,t') gives a section of this bundle. 

These geometrical constructions generalize to systems 
with n incoming and m outgoing channels, II. The first 
n rows of S span an n-dimensional complex subspace 
of C"+ m . The space of all such subspaces, called a 
Grassmannian, has a naturally defined 2-form, called the 
Kahlcr form p0| . Up to a constant factor, the Brouwer 
2-form equals the Kahler form. In addition, there is a 
canonically defined line bundle over this Grassmannian, 
and A equals the global angular form for this bundle. 

Statistics of weak pumping: Next we consider how 
a random scatterer transports charge when two param- 
eters are varied gently and cyclically. More precisely 
we consider the charge transported by moving along 
the circle X\ = ecos(r), Xi = esin(r) in parameter 
space. If e is small, then the charge transport is close 
to — Q(di, 82), evaluated at the origin, where dj are 
the tangent vectors associated with the parameters Xj. 
The vectors dj map to random vectors on U (2), which we 
assume to be Gaussian with covariance C . The problem 
is then to understand the possible values of the curvature 
ft applied to two random vectors. 

To do this, we first need to understand the statistics 
of 2-forms applied to pairs of random vectors, and to 
understand the geometry of the group U{2). 

Take two random vectors in R 2 , and see how much area 
they span. By random vectors we mean independent, 
identically distributed Gaussian random vectors whose 
components Xj have the covariance (XiXj) — CSij. The 
area A is distributed as a 2-sided exponential: 



sin 2 {d(f>) 2 + {d0f 



(14) 
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This is seen as follows. If the two vectors are X and 
Y, then the area is X1Y2 — X2Y1. X{Yi and —X2Y1 are 
independent random variables, and a calculation shows 
that their characteristic function is 1/ \Jk 2 C 2 + 1. Their 
sum is a random variable with characteristic function 
(k 2 C 2 + 1) j and so exponentially distributed. 

We parameterize the group U(2) by the angles 
(a, j,<p,6), as in (^J). A standard, bi-invariant metric 
on U(2) is 
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In this metric the vectors di are orthogonal but not or- 
thonormal. Unit tangent vectors are 
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■ d a , 
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d</,, eg = dg. (15) 



The volume form is sin(0) cos(0) da A d'y A d<f> A d0. The 
curvature 2-form, from Eq. (13), is 



Q = —2 sin(0) cos(0) d0 A (da + d(j>). 



(16) 



A scattering matrix is time reversal invariant if and 
only if t = t' . The space of time- reversal matrices is 
parameterized exactly as before, only now with (j) iden- 
tically zero. The volume form for the metric inherited 
from U(2) is cos(6)da A d'y A d0, and the curvature form 
is now f2 = —2 sin(0) cos(0)d0 A da. 

We are now prepared to compute the statistics of weak 
pumping, assuming that the S matrix is uniformly dis- 
tributed and that the tangent vectors to the space of S 
matrices are Gaussian random variables. This problem 
was studied by Brouwer Q , in the framework of random 
matrix theory, by which we mean that Brouwer posits an 
a-priori measure on the space of Hamiltonians. Random 
matrix theory is more powerful, in that the distribution 
of tangent vectors is fixed by the theory. The price one 
pays is that the analysis is also far from elementary and 
the results are, in part, only numerical. 

For systems without time reversal symmetry, random 
matrix theory posits that the S matrix is distributed on 
U(2) with a uniform measure. Since the conductance 
g is g cx \t\ 2 = sin 2 we have that dg cx sin 6 cos 0d0, 
proportional to the volume form: The conductance g is 
therefore uniformly distributed. 

A random tangent vector to U(2) is X — Xgeg + 
X a e a + X^e^ + X 1 e 1 , where Xj are Gaussians with 
(XjXk) = CSjk- The curvature associated with two ran- 
dom tangent vectors X, Y is, by Eq. ( p^ ) 



Sl(X, Y) = -2(x e W e - Y e Zg^j , 



(17) 
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Tr(dS* ® dS) = (d 7 ) 2 + cos 2 9 (da) 2 + 



where Wg — sin 6>Yq + cos0Y^ and Zg = sm0X a + 
cos9 X^. The variables Wg and Zg are independent, each 
with variance C. From Eq. (|l3|), the distribution of the 
curvature is exponential and independent of \t\. The joint 
distribution of curvature, w, and conductance, g ■■ 
is given by the probability density 

^e-^ 2C dcdg (18) 

with u> ranging from —00 to 00 and g from to 

For systems with time reversal symmetry, the S ma- 
trix is uniformly distributed on the t = t' submani- 
fold, with the metric inherited from U(2). The tangent 
vectors are now Gaussian random variables of the form 
X = Xgeg + X a e a + X-ye 7 , and the curvature is now 
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n(X,Y) = -2sm6(x e Y a -Y e X a y (19) 

Since the curvature depends on 8 the curvature and the 
conductance are correlated. The volume form indicates 
that ^Jg, and not g, is uniformly distributed. This fa- 
vors insulators. The joint distribution for curvature and 
conductance is 

1 e-M/tcjSnSfoj^ (20) 

This formula says that, statistically, good pumps tend 
to be good conductors; co/-Jg, rather than to itself, is 
independent of g. 

We have assumed, so far, that the variance C is a con- 
stant. There is no reason for this and it is natural to 
let C itself be a random variable. Given a probability 
distribution for the covariance, g?/i(C), one integrates the 
formulas (|l8|) and (|2(]) over C. One sees, by inspection, 
that in the absence (presence) of time reversal symmetry, 
lu (uj/^/g) is independent of g. Furthermore, the distri- 
bution of lu after integrating over g is smooth away from 
lu = 0, but has a discontinuity in derivative (log diver- 
gence) at uj = 0. In these qualitative features, our results 
agree with Brouwer's. However, the tails of the distribu- 
tion for large pumping may depend on the tail of d/i(C); 
While (|l^) and ( po|) have exponentially small tails, power 
law tails d[i(C) will lead to power law in the tails in lu. 
Since we do not determine dfi(C) we can not determine 
the tails. Using random matrix theory Brouwer deter- 
mined the power decay in u> || . 

Hard Pumping: Finally, we consider what happens 
for hard pumping. Here one can no longer evaluate the 
curvature at a point and multiply by the area. One needs 
to honestly integrate the curvature. Hard pumping was 
addressed by || who studied it in the context of ran- 
dom matrix theory and showed, using rather involved 
diagrammatic techniques, that pumping scales like the 
root of the perimeter. Here we shall describe a com- 
plementary, elementary result that holds provided the S 
matrix is a periodic function of the parameters. This is 
the case, for example, when the pumping is driven by two 
Aharonov-Bohm fluxes. 

With S(x, y) periodic in the driving parameters x and 
y, so is the curvature Sl(x,y) — ^mn e l ( mx + n v) _ Since 
the global angular form is also periodic, (Iqo = 0. 

The integral f\ x \ <R Q of the curvature on a large disc 
of radius R is, to leading order, 

V^V /2 sm (N(n,m)R-^), (21) 

where N(n, m) = \Jn 2 + m 2 . The charge transported in 
a cycle is proportional to the square root of the perime- 



ter (or the fourth root of the area) times a quasiperiodic 
function of R. This follows from the evaluation of the 
elementary integral J^ <R dx dy e^ nx+my \ which equals 
a Bessel function whose large- R asymptotic behavior is 
sin(Ai?-f ). FromEq. © one can determine the 

probability distribution for charge transport, (viewed as 
a random variable with uniform distribution on the ra- 
dius R). Eq. ( |2~l"| ) turns out to be closely related to a cel- 
ebrated problem in ergodic number theory: The Gauss 
circle problem [ p"l"| . 

This result does not directly apply to the pump stud- 
ied by Switkes 0, because the parameters they vary do 
not have built in periodicity. Nevertheless, it illustrates 
two features of pumps that have been observed experi- 
mentally. The first is that hard driving transports a lot 
of charge, with scaling that is sublinear, as in a random 
process, and the second that the directionality of hard 
pumping is essentially unpredictable. 
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